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MINIMIZATION OF THE TOTAL HEAT INPUT FOR MANNED VEHICLES 


ENTERING THE EARTH'S ATMOSPHERE AT HYPERBOLIC SPEEDS 

By Alvin Seiff and Michael E. Tauber 
Ames Research Center 


SUMMARY 


The use of conical entry bodies to control the radiative and minimize the 
total heat input to manned vehicles is investigated with an imposed accelera- 
tion limit on two types of hyperbolic shallow-angle entry into the Earth's 
atmosphere. The local Reynolds number is limited to retain laminar boundary 
layer. The equations permit any desired values of the acceleration and local 
Reynolds number limits and of the body radius , body density, and heat -shield 
material properties. Numerical examples are given for a range of Reynolds 
number limits, three body radii, and two heat-shield materials. 

The dependences of the convective and radiative heating rates on air 
density, velocity, cone angle, and base radius are written and integrated over 
the trajectory in an expression which is particularly simple for nearly hori- 
zontal flight. The major reduction in convective heating due to transpiration 
from the heat shield is included, although somewhat uncertainly for the higher 
speeds. The heat inputs are expressed as fractions of the vehicle kinetic 
energy at entry to obtain a dimensionless figure of merit for comparing cases 
of different vehicle mass or entry velocity. Trajectories analyzed include 
constant altitude deceleration and constant Reynolds number descent. The 
latter is the trajectory of minimum convective heat input for a given Reynolds 
number limit, and of minimum total heat input in cases where the radiative 
contribution is small. Aerodynamic lift required to follow these trajectories 
is discussed. 

For any specified set of entry conditions, the cone angle has an optimum 
value which diminishes with increasing entry velocity. The optimum cones gen- 
erally have small total radiative heat input compared to convective. For cone 
angles larger than optimum, the radiative and total heat inputs can become 
many times greater than the optimum. The total heat input energy fraction is 
a sensitive function of body size on acceleration-limited entries. It is 
smallest for large bodies of the order of 3- to 4-meters base radius - for 
example, less than 1 percent of kinetic energy at entry for velocities up to 
about 24 km/sec. The corresponding mass loss may be less than 20 percent of 
total mass at entry. With small bodies or very low Reynolds number limits, 
the optimum energy fraction is increased several fold. While the heating may 
be manageable, the requirements for lift-drag ratio and acceleration limiting 
appear to become difficult to satisfy aerodynamically with optimum cones on a 
single pass at entry velocities above about 25 km/sec. 



INTRODUCT ION 


As space flight missions go to progressively greater distances from the 
earth, the speed of entry into the atmosphere on return to earth tends to 
increase. Thus, intercontinental ‘ballistic vehicles and near-earth satellites 
enter the atmosphere at 6 to 8 km/s; returning lunar vehicles, at 11 km/s; 
and vehicles returning from the near planets, at speeds as high as 20 km/s. 
Still higher entry speeds are likely to arise in connection with future 
missions not yet considered in detail. 

At the lower speeds, the entry heating has been readily controlled by the 
use of configurations which direct most of the vehicle kinetic energy into 
heating the atmosphere (ref. l), and by the use of ablation shields which 
block and absorb the energy converted to heat in the boundary layer. However, 
the total energy to be dissipated (per unit of vehicle mass) increases with 
the square of entry velocity, and, as will be shown, the fraction of the 
vehicle kinetic energy deposited as heat in the body also tends to increase 
with entry velocity. Successful solution of the heating problem at high entry 
velocities depends on maintaining this heating energy fraction small. 

The principal mechanism for body heating at velocities less than 10 km/s 
is boundary-layer convection. At higher velocities, the gas in the shock layer 
may attain temperatures well in excess of 10,000° K and radiate substantial 
amounts of energy at optical and ultraviolet wavelengths (ref. 2). This radi- 
ation falling on the body tends to become the predominant heating mechanism 
for blunt bodies and increases the fraction of the vehicle energy converted to 
body heat input . 

While the intensity of radiative heating at any wavelength cannot exceed 
that of a black body at the temperature of the shock layer and, for some low 
altitude flight conditions, may approach that limit, such heating rates would 
generally be considered unacceptably large. For example, at the modest shock- 
layer temperature of 65OO 0 K, the black body intensity is 10 5 kW/m 2 . Heating 
rates are also limited by the energy available in the incident flow. For the 
portion of the trajectory where the deceleration is large, the available flow 
energy is also large ((l/2)pV 3 A = ma x v/Cp), and substantial fractions of it 
can be radiated toward the body. Thus, these two limiting effects do not 
prevent the radiative heating from assuming very large values. 

It has been shown in references 3 through 6 that radiative intensity is 
reduced for bodies, such as cones, which have oblique bow shock waves. As is 
well known, the oblique shock wave is equivalent to a normal shock wave at the 
lower flight velocity U = V sin 0 W . Furthermore, the radiative emission 
behind a normal shock wave varies as a high power of the velocity (refs. 2 and 
5) so that sweepback of the shock wave offers a powerful means of reducing 
radiative input. Bodies having oblique shock waves, however, exhibit smaller 
drag coefficients than blunt bodies, and this results in increased convective 
heat input, as will be shown. Therefore, an optimum cone angle will be found 
which minimizes the sum of convective and radiative inputs. The analysis of 
reference 5 defines these optima for ballistic entry. It was found that the 
fraction of entry kinetic energy that enters the body as heat may be kept 


2 



below 0.01 at optimum for entry velocities to 30 km/s, if the trajectory is 
limited to a peak Reynolds number below the transition Reynolds number. It 
follows that such high-speed entries might be made with a relatively small 
heat-shield mass loss. 

The vehicle accelerations on ballistic entry at such velocities are large, 
hundreds to thousands of times the earth 1 s gravitational acceleration, and are 
far in excess of human tolerance. The question of optimum bodies for manned 
entry at hyperbolic speeds therefore remains to be treated and is the present 
subject. The restriction to suitably low levels of acceleration is obtained 
by choice of trajectories of shallow angle, sometimes called grazing trajec- 
tories. Two kinds of grazing trajectories are treated - constant altitude 
decelerations and constant Reynolds number descents. The limitation imposed 
on entry by the lift required at supercircular speeds is considered. Optimum 
cone angles are defined, and the fractions of entry kinetic energy that go 
into heating the body are calculated. Simplifying assumptions are made, so as 
to obtain generality, and a range of entry parameters is considered. The 
heating analysis follows that of reference 5* The problem of ablating the tip 
during entry is analyzed to determine if an effectively sharp conical body can 
be maintained . 

Preliminary results of the analysis have been given in references 8, 
and 9 • 


NOTATION 

a entry vehicle acceleration, m/s^ 

a shock -wave area, m 2 

A body frontal area, m 2 

B ballistic parameter, Cpp 0 A/3m sin /p; 

Cp gas specific heat at constant pressure, j/kg °K 

C c laminar convective constant (eq. ( 17 )), (s/m) J m° #5 

Cp drag coefficient, dimensionless 

C e equilibrium radiative constant (eq. (10)), (w/m 3 ) (s/m) ^ 

Cjj heat- transfer coefficient, dimensionless 

C n nonequilibrium radiative constant (eq. ( 13 )), (w/m 2 ) (s/m) S 

D drag force, N 


3 



power radiated from gas in shock layer , W 

convective heating correlation constant (eq . (Al)), (m/s) 0,36 
acceleration due to gravity, m/s 2 
local static enthalpy, J/kg 
streamline total enthalpy, j/kg 

heat input integrated over body and over time, J 
defined in equation ( 13 ) 

laminar convective heating integral (eq. ( 29 )), (m/s) J+2 
velocity exponent (eq. (l6)) 
constant (eq. (35)) 

pV, a constant on constant Reynolds number trajectories, kg/m 2 s 

a /t, (s/m ) 2 

lift force, N 

entry vehicle mass, kg 

Nusselt number (eq. (A2)) 

density exponent (eq. (10)) 

cone surface pressure, N/m 2 

velocity exponent (eq. (10)) 

local laminar convective heating rate (eq. (A2)), w/m 2 

free-stream dynamic pressure, appendix C, N/m 2 

body base radius, m 

body nose radius, m 

gas constant, j/kg °K 

Reynolds number based on base radius and free-stream air properties 

Reynolds number based on slant length and properties at boundary- 
layer edge 



Re ^ maximum allowable value of Re G 

Rp radius of planet, m 

s velocity exponent (eq. (13)) 

5 slant length of cone 

t time, s 

T temperature, °K 

u air velocity in the shock layer, relative to the body, m/s 

U component of flight velocity normal to shock wave, m/s 

Ui ^ 2 velocity of change in the equilibrium radiative constants, 13^700 m/s 
v volume, m 3 

V flight velocity, m/s 

V s near-earth satellite velocity, 7900 m/s 

W entry vehicle weight, N 

x range measured along flight path, m 

x n rate of axial regression of the tip due to ablation, m/s 

y alt itude , m 

Z moles per original mole 

a transpiration effectiveness parameter 

Of- trim angle of attack 

P reciprocal of atmospheric scale height (eq. (Cl 8 )), m ” 1 

7 flight-path angle below horizontal 

6 shock-wave standoff distance (eq. (e 4 )), m 

£ energy to ablate unit mass of heat shield, j/kg 

H 

r) — — , dimens i onles s 

mV E /2 

G c cone half -angle 


5 


conical bow wave half -angle 


M- 
P 
P 

pb 

a 

+ 

C) 

Subscripts 

c laminar convective 

cl collision limiting 

c Q laminar convective in absence of ablation 

e equilibrium radiative 

E at entry into the atmosphere 

1 limiting or maximum permissible value 

min minimum 

n nonequilibrium radiative 

o earth sea-level condition 

opt opt imum 

st stagnation point 

w evaluated at vehicle surface 

x horizontal component, or component along the flight path 

y vertical component, or component normal to flight path 

e evaluated at boundary-layer edge 

2 shock-layer properties 

Stream properties without subscript are evaluated in the free stream. 
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gas viscosity, kg/m*s 
air density, kg/m 3 
p/po 

body density, kg/m 3 

value of ^ for V °o ( e q. (l8)) 

V^co 

derivative with respect to time 



ANALYSIS 


Heating Equations 

The rate of heat energy input to a vehicle flying in the atmosphere may 
be expressed as a fraction Ch of the kinetic energy pV 3 A/2 of the air 
intercepting the body per unit time 

(i) 

To satisfy the conservation of energy, the heat -transfer coefficient Ch 
defined by equation (l) cannot exceed 1.0. The work converted to heat by the 
total drag force per unit time 


DV = Cd 


( 2 ) 


results partly in heating the atmosphere and partly in heating the entry 
vehicle. That is, dH/dt < DV, so that Cn/Cp is less than 1 and represents 
the fraction of the work done by the drag force which goes into heating the 
vehicle . 

The total heat input on a trajectory is obtained by integrating 
equation (l) 



Equation ( 3 ) is evaluated in reference 5 for rectilinear ballistic entry into 
an exponential atmosphere. For other trajectories, it can be integrated if 
the dependence of Ch, P, and V on t is established. Since Cg; can be 
expressed as a complicated function of p, V, and for a given body shape, 
relations between p, V, and t which permit the expression of the integral in 
terms of one variable are sought. These relations are provided by the 
trajectory equations. 

The rate of velocity change along the flight path may be written 



W 

Sketch (a) 


~ m It = Cd ~ “s sin y (*o 

For trajectories within the atmosphere, with the 
drag force comparable to or larger than the weight, 
and at small values of 7, this becomes, very 
closely. 
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from which 


-m f = ° D 


pV 2 A 

2 


at = - -S_2 dv 

CdA P v 2 

Ulus, for flight at or near zero path angle 7 , equation (3) becomes 


H = 



Cg 

-S mV av 
CD 


( 5 ) 


(6) 


•which depends on p only through the dependence of Cg on p. For reference, 
it is informative to consider the case Cg/Cp = constant. Then, with constant 
vehicle mass, equation (6) may be integrated to obtain H = (Cg/Cp) (mVE 2 /2) . 
Here, Cg/Cp is seen to be the fraction of kinetic energy at entry which is 
converted to body heat. With Cg variable, the total heat input H may still 
be expressed as a fraction g of the kinetic energy of the vehicle at entry 

h ( 7) 

2 

where g is a weighted mean value of Cg/Cp, defined by combining equations (6) 
and (7) to obtain, for a constant vehicle mass. 


g = 



£h _v_ _v_ 

Cp Ve Ve 


(8) 


It is seen that g may be made small by minimizing the area under the curve 
(Ch/C D )V as a function of V. Note that equations (7) and (8) show that the 
total heat input is inversely proportional to drag coefficient. 


The mass lost by ablation may be written, for cases where the heat stored 
is small compared to that resulting in ablation. 

An _ ^E 2 

m E 2 £ 

where £ is the energy absorbed per unit mass of material ablated. 


(9) 


The Heat -Transfer Coefficient 

The heat -transfer coefficient Ch is a function of the flight velocity, 
stream density, body shape and size, atmospheric composition, and nature of 
the boundary layer. Even when the boundary layer is laminar, the heat trans- 
ferred by convection at the speeds considered is very large; while with a 
turbulent boundary layer, the convective heat input integrated over the tra- 
jectory is from 5 to 10 times larger (ref. 5 )* The present analysis seeks to 
define conditions of minimum total heat input, and these will clearly corre- 
spond to those which maintain a laminar boundary layer. It is assumed that 
laminar flow can be maintained by restricting the local Reynolds number below 
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some specified limit. The appropriate critical value or values for this limit 
are presently unknown for conditions of high total enthalpy flow over an 
ablating surface. Therefore , the value to be chosen is treated as a parameter 
of the analysis, and the analysis is restricted to a consideration of the 
laminar flow boundary layers. 


The heat-transfer coefficient is necessarily treated in three parts, the 
laminar convective contribution, Ch C j the equilibrium radiative contribution, 
CHe> and the nonequilibrium radiative contribution, CHn* The latter two 
are assumed to be given by equations from reference reproduced below. 


The radiative power emitted per unit of shock-layer volume at equilibrium 


is 


= CeU^ 
dv 


( 10 ) 


where U = V sin 0-w is the component of velocity normal to the shock wave and 
p = p/p Q . Values of the constant and exponents are , for air , 

, , / NlS^S^ 

r = 6.l4 _W_ [ £ ) 
ei 10 49 m 3 W 


and 


q. x = 15*^5 


p = 1.80 


> U < 13,700 m/s 


-e 2 


= 


6.4-4 jw_ 

]_0 6 m3 Vm. 

5-05 


5 .05 


/ U > 13,700 m/s 


(Ha) 


(lib) 


with p 0 
written 


p = 1.80 

1.225 kg/m 3 . The equilibrium radiative component of C g; may be 


„ _ TT (q.-3)_(p-i)ru tan 2 0 W - tan 2 0 C , . a 

^ ^ P £ 3 £ n3 Uln ^ 


Po 


( 12) 1 


Values of the wave angle for cones in high velocity flow are given in 

figures 35 and 36 of reference 5j based on the analysis in appendix B of that 
reference . 


^or considering the dimensional consistency of equation ( 12 ), note that 
in the mks system the W/m 3 is identical to the mechanical unit kg/(m)s 3 . 
Similarly, in equations ( 13 ), (1*0, and (15) the w/ m 2 is identical to the kg/s 3 . 
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The nonequilibrium. radiative power unit of shock-wave surface is assumed, 
following arguments given in reference 2, for example, to be a function of the 
velocity U only, above a certain ambient density, 


where 


and we have used p , 

cl 


dE = 
da 


C n U' 


£_Y 

Pc lJ 


G n 


1.15 X 
]_qH m2 



10 “ 3 as 


in reference 5* 


s = 4.1 


(13) 

(14) 


In the "collision limited" 
the nonequilibrium component of 
from this is 


c H n “ G n 


V s-3 

P 


regime i = 1, while for p > P c ^, i = 0, and 
the heat -transfer coefficient which follows 


tan 0 W \ 2 / . „ . 

7 t 2 -) (sin e w ) 

tan 6 C / 


s-i 



(15) 


The convective heat-transfer coefficient may be written, based on a recent 
analysis of equilibrium real air boundary-layer heat transfer to pointed cones 
( ref . 10 ) , 


with 


Ch Co 1 — ^ sin 2e c 

*7^1 bp 

( 16 ) 

/ s 1 /e 

n 0.84 f s \ 0.5 

- 10 4 y - 

(17) 


J = l / 6 


It will be noted that this coefficient and exponent differ from those given in 
reference 5- Appendix A derives equations (l6) and (17) from the correlation 
given in reference 10, and compares the values of Ch Co "with those of 
reference 5. 2 


The convective heat transfer in the presence of ablation is reduced by the 
transpiration cooling effect of the ablation vapors. Following reference 5 we 
assume 


t = 


C H C 

C H Co 


1 - p 
1 + KcV 2 


+ P 


(18) 


Values of C c and j derived in appendix A for T w = 3OOO 0 K, and 
between 0.1 and 1.0 atmosphere are I.l6xi0 -4 and 0.l4, respectively. For 
velocities between 10 and 20 krn/s, values of Ch Co given by these constants 
and those of equation (17) differ by less than 5 percent. The choice of 
j = l/6 is advantageous because it permits integration of the convective heat- 
ing over the trajectory in closed form, as will be shown. 
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where Kc = a /£ , a, is a laminar transpiration effectiveness parameter that 
depends on the molecular weight of the ablation products (ref. ll), £ is the 
energy per unit mass (mechanical energy units) required to heat, dissociate, 
and vaporize the ablation material, and o is an assumed asymptotic value of 
Ch c / c Hc 0 for velocity. 

With a = 0 , equation (l8) gives t = 0 at infinite velocity. This cor- 
responds to infinite effective heat of ablation (a concept not used herein) at 
infinite total enthalpy. Provision of the term a as an asymptotic limit of 
for V -+ oo is based on experimental results from Ames Research Center arc 
jet facilities (ref. 12) which suggest that the factor ^ does not go to zero 
but to a small fractional value of the order of 0.1. Since the behavior of 
ablation shields at the high velocities considered here is not well estab- 
lished, the effect of varying the asymptote a from 0 to 0.1 will be shown. 


It may also be remarked that for charring organic ablators, the true heat 
of ablation £ is not a constant, being a function of the degree of dissocia- 
tion of the gaseous products and, therefore, of the temperature of the ablat- 
ing surface. Hence, it depends on the heat-transfer rate and distribution 
between the convective and radiative contributions. Since the study of such 
materials is quite specific to the material chosen and too complicated for 
parametric treatment, it was not attempted here. The materials chosen for the 
examples were a low temperature organic ablator. Teflon, and a material with 
a large heat of vaporization, quartz, which is assumed to vaporize and not run 
off the surface as a liquid, (in practice, quartz does suffer liquid runoff, 
particularly at moderate heating rates.) The choice of these materials is 
intended to be illustrative, and does not imply a judgment by the authors 
that they represent leading contending materials for high speed entry. A 
third material, graphite, was included because it appeared to offer advantages 
for protection of the conical tip. The assumed properties of these materials 
are given in the following table. 



a 

(m/s) 2 

, Kc, 

(s/m) 2 

Teflon 

0.26 

2 .2X10® 

12.0X10 -8 

Quartz 

0.24 

16 .0X10® 

1.5X10" 8 

Graphite 

0.42 

67 .0X10 6 

0 .63x10” 8 


For comparison, the heat of ablation £ of the charring ablator, phenolic 
nylon, at a high surface temperature, 3500 ° K, has been estimated to be as 
high as 28 X 10 6 (m/s) 2 (ref. 13) . Thus, its mass losses, exclusive of chemical 
and mechanical erosion, should be somewhat smaller than those of vaporizing 
quartz as presented herein. For graphite, a is again chosen to follow the 
molecular weight dependence suggested in reference 11, and the heat of abla- 
tion is taken from data in reference l 4 . 
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Integration of the Heat-Transfer Equation 


In terms of the components of heat input, Ch = Ch c + OHe + CHru equa- 
tion (8) becomes 

PiXji (19) 

' J 0 C D V E V E J 0 C D V E V E J 0 C D V E V E 


or 


T) = Tic + % + Tin 


(20) 


From (19), (20), and (12), with Cp = 2 sin 2 0 C , corresponding to the Newtonian 
drag coefficient of a cone at zero angle of attack, 3 

ne = fV-’p 1 " 1 a tana f 9 ° (sln ^ f 4 f (21) 

Po 3 tan 3 0f> sin 2 Qr. ^E V E 


Assuming that the cone angle and base radius are constant during entry and 
that the bow wave angle is independent of speed, we write, noting that 

v/v E = u/u E , 


'He 


r^ tan 2 0 W - tan 2 0 C 
Po P 3 tan 3 0 C sin 2 9 C 


. -3 


C ei V E 1 (sin 0 W ) 1 


E ^(±) a 
Vu E / 


u 

Ue. 


+ Ce 2 v E 2 (sin 9 v ) q ~ 2 f df U 


Ul,2 

U E 


U: 


Ue. 


(22) 


If a constant altitude is maintained during entry, p is constant and 
equation (22) may be integrated to obtain 


Tie = *b 


where 


pP -1 (tan 2 0y - tan 2 0c) sin 0y 
Po^Vv 2 3 tan 3 0c sin 2 0c 


C~ U‘ 




C e U 




'I 1 »2 + ^2 1>2 


q x - 1 q 2 - lL VUi,a 


U E 


q*- 1 


- 1 

(23) 


and 


q. -1 
C U 1 

. , e 1 A; A-. = 2.55X10 
9i “ 1 

PL0-1 


10 _J_ = kg 

m 4 (m-s) 2 


C e U 

e 3 

q 2 - 1 


e 2-Jo 2 _ = 9. 02X10 10 

m 4 


' ^The question of how such a cone can develop the lift force generally 
necessary for hyperbolic grazing entry will be discussed in a later section. 


12 


I 



Alternatively, if the Reynolds number is maintained constant along the 
trajectory, pV is approximately constant and 


Tie 


V E ‘ 


Po 


where 


and 


> w - tan 2 0 c )(sin 

q, ~p 
C e U^/ 

. -1 

tan 3 0c sin 2 0c 

q x - p 


q 2 -P 


+ - 


q P - p 


u E 

JJi,a 


q 2 -p 


-i 


(24) 


q., -P 

-1 

<Li ~ P 
.cu-P 


32X10 7 4 


m' 


1.80 


Ce 2 U l ? 2 , _ 5.5IXIO 7 JL (1 

q - p m3 V» 


1.80 


In equations (23) and (24), when U E < U^s, Ui 2 is to he replaced by Ug. 

The nonequilibrium contribution to t) is written by use of equa- 
tions (19)^ (20), and (15) 


Tin = c nV E " 


s-3/tan 0 V 
£an 0 Ci 


(sin e w ) s-1 r 1 (v/v E ) s " 2 
sin 2 0 C vi 0 P 



which, for constant altitude entry, becomes 


— ©w\ 

nn ” p(s - 1) vtan 0 C / 


(sin ey) 83 " 1 / pY 
sin 2 0 C \P c is 


(25) 


(26) 


and for constant Reynolds number entry 


nn = CnVB*" 


s _ 3 /tan 0. 
ban Q r 


w 


(sin 0 V ) S ~ X 
sin 2 0 C 


PczPo^ 3 - 1 ) 




V E 

( pV)s 



(2Ta) 

For the cases of interest here, it is generally permissible to simplify equa- 
tion (27a) to 


n _ °nV E 

’“'TFi: 


S-2 


tan 0. 
±an Or 


w 


(sin 0 V ) S ~ X 
sin 2 0 C 


(27b) 


where pV is a constant of the motion. 

The laminar convective portion of 7 is written by use of equations (l6), 
(18), (19), and (20) 
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dV 


(28) 


= 


C c ^sin 20c 
Vg 2 >Jrb sin 2 6, 


C 4> NL + KqV 2 / VW 


This integral cannot be evaluated in closed form for arbitrary noninteger 
values of j. However, it can be readily evaluated numerically. For constant 
altitude entries, the constant s/W is taken outside the integral; the integral 
is then a function of Ve and the ablation material constants K c and 0 only. 


I c (V E ,Kc,0) 


= f* (l-a + „ j 
1 Vl + K C V 2 J 


yj +1 dV 


( 29 ) 


For given material properties,, it need be evaluated one time only, as a func- 
tion of Ve- Values of this integral normalized on Ve j+2 are plotted in 
figure 1 for the material properties previously tabulated for Teflon and vapor- 
izing quartz and for o - 0 and a = 0.1. 
It can be noted that the choice of 
value for cr is not a critical one for 
vaporizing quartz within the velocity 
range considered, but is very important 
for Teflon at the higher velocities. 



Quartz 


Teflon 


On constant Reynolds number tra- 
jectories, the integral can be obtained 
in closed form for the special case 
j = l/6 by the method shown in 
appendix B. 


Heating Below Satellite Speed 


Figure 1. - The laminar convective heating 
integral for various assumed, ablator 
properties. 


Below satellite speed, only con- 
vective heating is important . Since 
this speed range contributes only a 
modest fraction of the total heating during the entry, no effort was made to 
optimize this part of the flight path. (Extensive work dealing with the prob- 
lem of entry at satellite speed has been reported in the literature.) 


For constant Reynolds number trajectories it was assumed that the 
constant Reynolds number flight path was maintained to terminal velocity. In 
actuality, as will be shown, such a trajectory can be followed only to the 
point where Vy = V. However, this occurs at very low speed - a few hundred 
feet per second for the larger bodies, a few thousand feet per second for the 
smaller ones - and the heat input below this speed is an insignificant part 
of the total, so no appreciable error results from the assumption. 

For the constant altitude trajectories the constant altitude flight paths 
were followed to somewhat below satellite speed where they were matched in 
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velocity and altitude to zero-lift ballistic decay trajectories originating at 
satellite speed above the atmosphere . A small discontinuity in 7 that occurs 
at the matching point was ignored. The time history of the coordinates of this 
kind of flight path are tabulated in reference 15 * The convective heat- 
transfer relation, equation (28), was applied to the appropriate portion of 
these trajectories. The maximum acceleration is 8.3 g> making the satellite 
speed descent trajectory suitable for manned vehicles. 


Trajectory Analysis 

The choice of trajectories was important to the analysis. The constant 
altitude deceleration and constant Reynolds number descent trajectories were 
selected to simplify the heating analysis as well as to permit studies of 
parametric variations in vehicle size, density, cone angle, etc., on similar 
flight paths. At the same time they are representative of the heating on 
shallow angle entry. Trajectory types considered in some earlier work, such 
as constant lift-drag ratio trajectories, are, by comparison, complicated by 
altitude excursions (skipping motions), lack generic similarity, and can be 
analyzed only numerically. One of the types chosen here, the constant Reynolds 
number descent, is thought to present the minimum heat input to bodies of 
optimum cone angle. As such, it is a useful standard of comparison. 

The objectives of the trajectory analysis were to determine the feasibil- 
ity of the trajectories for real applications, to supply needed constants for 
the heating analysis (free-stream density, primarily), and to obtain closed- 
form expressions for the usual trajectory variables of range, acceleration, 
velocity, etc., as functions of time. The feasibility of the trajectory is 
largely determined by the lift -drag ratio required to follow it without exceed- 
ing the constraints on acceleration and Reynolds number . The constraints 
amount to a specification of minimum permissible altitude, in the case of the 
constant altitude deceleration, and of the altitude at the point of entry onto 
the constant Reynolds number descent, in the other case. 

It was assumed that the lift-drag ratio would be generated by a geometric 
modification of the cone, such as those described in appendix C. Roll oscilla- 
tion at the trim attitude provides modulation of the vertical component of 
L/D, which is required, in general, to be a decreasing function of time. 

The principal trajectory equations are given below and developed in 
appendix D . 

C onstant altitude deceleration .- The lift-drag ratio required is at every 
instant a function of the total acceleration and the flight velocity 



where the total acceleration, the vector sum of that along the path and across 
the path, is given by 


15 


1/2 


a 

g 



PV. 


2 “I 
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(31) 


and w/CdA may be written in terms of the cone angle, mean body density, and 
base radius . 


W = Pb^bg 

Q[)A 6 sin 2 0 q tan 0 C 


(32) 


When the total acceleration is limited to a specified value, a i, the 
corresponding L/D is given by equation ( 30 ), and the ambient density at the 
flight altitude may be obtained from equations ( 31 ) and ( 32 ) as 


P 


Pb r bg 

3V S 2 sin 2 0 C tan 0 C 


( a 7,/g) 2 

-(Ve/Vs) 4 


li/ 2 


1 


(33) 


If a Reynolds number limit rather than the acceleration limit determines the 
altitude , the density is given by 


where 


P = Re 2 


tan Q r 


k V E r b 


Re e 

k ^ — tan 0 C 


(3*0 

(35) 


is approximately constant and equal to 2 for U < V s (appendix D). For 
U > V s , k takes on smaller values, going to about 1»3* 


The velocity-distance, velocity- time, and time -distance relations are 
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( 36 ) 


and 


V = 


V E 


1 CpA 

1 + 2 ~m~ P v Et 


1 Cr,A CtsA 

e 2 m px = i + l pVjrb 

2 m 


( 37 ) 


(38) 


respectively, where x is arc distance along the orbital path. Fixing the 
maximum acceleration and the entry velocity can be shown to determine (CpA/m)p. 
Hence , for given maximum acceleration and entry velocity, the x - t and V - t 
relations are universal. 
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Constant Reynolds number descent .- Appendix D shows that to maintain a 
constant Reynolds number requires a constant rate of descent, Vy = constant, 
given by 

v _ PdA- (pv) 

Vy " ~ ~2f~ 



where pV is a constant and (3 is the reciprocal of the atmospheric scale 
height. The angle of descent is not constant, however, being defined by 


Vy 

Bin 7 = y- 




where V decreases with increasing time . 4 The angle of descent increases 
until, in some cases, Vy = V and 7 = 90°. At this point, the motion is 
governed by the equations of simple ballistic descent and it is no longer 
possible to maintain the Reynolds number constant. 

The lift-drag ratio is governed by equation (D3^), which may be simpli- 
fied, for the cases of most practical interest, by use of small angle approx- 
imations for the sine and cosine of the flight-path angle to obtain 


L = 1 - (V S /V ) 2 + Vy 

D PRp( Vy/v ) v 


(M) 


where Rp is the radius of the earth or other planet whose atmosphere is 
entered. The total acceleration may be written 


a = ( V_ 
g \y s 


(pv)v 


S 


:] * K® * >]7 


2 j 1/2 


-2(V/V S )(W/CDA)_ 

which may also be put in terms of the rate of descent to obtain 

\ 4 _ /„ T7- \ 2 / \ 4_ i 1/ 2 


I = [(%) ( pRp)2 + ® [( pr p )2+ 2 pr p ] + (£) _ 


(42a) 


(42b) 


Finally, the relations between velocity and time, velocity and distance, and 
time and distance are, respectively. 


4 The heating occurs predominantly at velocities greater than Vg /k and, 
in this range, sin 7 is smaller than 4 sin 7 e (eq. (40)). The magnitude of 
7 e is defined by equations ( 39 ) and (40), and it is found that Vy/Vs is, 
for cases of interest, a very small number, corresponding to angles smaller 
than a degree. Hence, the use of small flight-path angle approximations for 
L/D is justified, and the neglect of the weight force component in equation (4) 
is also valid. 
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where, again, pV 


in equations (^3) to (^5) is a constant. 


PRESENTATION OF RESULTS 


( 43 ) 

( 44 ) 

( 45 ) 


Trajectories 


Constant altitude .- Ambient densities and flight altitudes dictated by 
the acceleration and Reynolds number limits are shown in figure 2 as a function 
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of entry velocity for representative 
values of the governing parameters. The 
- =13,500 kg/m 2 steeper family of curves is the 

acceleration-limited set. The inter- 
sections of this family with the 
Reynolds number limited family define 
points of transition from Reynolds num- 
ber to acceleration governed altitudes. 
Consider the two dashed curves, for 
example . At entry velocities below 
1.66 V s , the Reynolds number requires the 
lower density and governs, while above 
this velocity, the acceleration limit 
governs . For a 3O 0 cone with a base 
radius of 3*05 m, the Reynolds number limits corresponding to the three curves 
are 10 7 , 5^10 6 , and 10 6 , respectively. For the same radius and cone angle, the 
acceleration curves represent body densities of 6^0, 320, and 160 kg/m 3 , top to 
bottom. For these conditions, then, a Reynolds number limit of 10 6 would 
govern altitude selection for the entire range of entry velocities considered. 



Figure 2.- Ambient densities and altitudes 
dictated "by acceleration and Reynolds 
number limits. 


Entry may take place at higher altitudes than those calculated, but not 
at lower altitudes without violating the imposed limits on acceleration and 
Reynolds number. The altitudes shown, therefore, represent undershoot limits. 
The overshoot limits, not shown, depend on the available l/d. The undershoot 
boundary, however, is the altitude of interest for this analysis, since it 
minimizes the total heat input for optimum bodies. 5 


5 Guidance limitations make it necessary that real vehicles be able to 

withstand entry over a finite corridor. They must therefore be designed to 

accept somewhat larger heat inputs than the minimum considered here. 
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The lift-drag ratios needed for maintaining constant altitude flight are 
shown in figure 3 as a function of velocity for several levels of total accel- 
eration. In general, the maximum L/d is required at the highest velocity, 



that is, at entry in our analysis. A lift-drag ratio of 1 is sufficient to 
permit a 10 g limit constant altitude entry at a velocity of 2.76 V s . At 
higher entry velocities, the L/D required rises rapidly, or, alternatively, 
the acceleration limit must be increased. For any acceleration lev el, the 
L/D curve is asymptotic to a vertical line at a velocity V;g/V s = ^/a/g. This 
is the velocity for which the centripetal acceleration is the total accelera- 
tion. The entry body is then flying above the atmosphere, but still requires 
lift to hold it in a circular path; hence, L/d must be infinite. 

When the altitude is governed by a Reynolds number limit, or is higher 
than the minimum altitude permitted, the acceleration is smaller than the 
acceleration limit. Figure 3 is then entered at the appropriate level of 
acceleration to define L/d. A consequence of such a situation is that high 
lift-drag ratios may be required at the lower entry velocities. 
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The variations of velocity, range, lift-drag ratio, and acceleration with 
time after entry are shown in figure k for an entry velocity of 2.5 V s and an 



t.min 

Figure 4.- Representative characteristics of an acceleration-limited constant altitude 

trajectory, a z /g = 10, V E /V S =2.5- 


acceleration limit of 10 g, attained at entry. The curves are terminated at 
satellite velocity. The range covered to this point is 4600 km, and the 



0 10 20 30 40 50xl0 3 

^b r b’ kq/m 

Figure Descent rates for constant 

Reynolds number as a function of body 
density and size# 


elapsed time is roughly 6 minutes. The 
lift-drag ratio and total acceleration 
decrease continuously with time, the 
former going to zero at satellite 
velocity. 

The curves in figures 2 and 4 are 
examples. Other cases of interest may 
be readily obtained from the trajectory 
equations . 

Constant Reynolds numb er. - The 
descent rates characteristic of con- 
stant Reynolds number trajectories are 
plotted in figure 5 as a function of 
the body density times the base radius 
squared, for two cone angles and two 
values of the Reynolds number. The 
equation given in the figure is an 
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alternate form of equation ( 39 ) which incorporates equations ( 32 ) and ( 35 )* 
Large or dense bodies tend to have quite low rates of descent to maintain con- 
stant Reynolds number. The assumptions cited earlier concerning small angles 
of descent are amply satisfied by the conditions in this figure. 

The lift -drag ratios required are plotted in figure 6 as a function of 
flight velocity for various rates of descent. At Vy less than about 25 m/s, 
the required lift-drag ratio starts to become large, and may exceed that which 
is available. The accelerations, on the other hand, are minimized by use of 
small descent rates (fig. 7 ) and 25 m/s is, fortuitously, again near the bound- 
ary for acceptable levels of acceleration at the higher entry velocities. 



v s 


Figure 6.- Lift-drag ratios required on 
constant Reynolds number descents. 



v s 

Figure 7»- Acceleration on constant 
Reynolds number descents. 


Thus, to satisfy the requirements on 
both L/D and acceleration, values of 
Vy must lie in the vicinity of 25 m/s. 
A band in figure 5 indicates this 
acceptable region. Values of body den- 
sity and base radius to which it corre- 
sponds depend on the cone angle and 
Reynolds number chosen and may be calcu- 
lated from the equation in this figure. 

An example of a constant Reynolds 
number descent that satisfied reasonable 
requirements on acceleration and L/d 
is shown in figure 8 . The body size 
and density give a value of the param- 
eter p^r -^ 2 of 2980 kg/m. For a cone 
angle of 30° and a Reynolds number of 
5 million, the descent rate, according 
to figure 5> is 27 m/s. The trajectory 
variables in figure 8 show a descent to 
satellite velocity in 4 minutes from an 
entry velocity of 2.5 V s with a maxi- 
mum acceleration of 9*8 g, a maximum 
L/D of 0.7, and a flight range of 
3150 km. The angle of descent is 0.08° 
initially, increasing to 0 . 2 ° at satel- 
lite velocity. Following a trajectory 
with such precision may, in practice, 
present a problem in flight-path con- 
trol. In principle, however, such 
trajectories appear feasible. 


Heating 

Constant altitude entries . - The 
total heat Input energy fraction to 
cones of large base radius with 30 ° a nd 
35 ° half-angles is shown in figures 9 as 
a function of entry velocity. For the 
conditions cited, the total heat input 
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ranges from 0.4 to 2.5 percent of kinetic energy at entry, increasing with 
increasing Ve^ and depending also on the cone angle and base radius. 

Laminar convective heating predominates in the case of the 30° cone with 
a base radius of 3 . 05 m (fig. 9( a ))- The equilibrium radiative heating is 
kept small by the obliquity of the shock wave. The nonequilibrium radiative 
contribution is essentially negligible. 

For the 35° cone with a base radius of 4.57 m (tig. 9(~b)) , however, 
equilibrium radiative heating becomes important at the higher entry velocities 
and, at V E = 3 V s , contributes more heat input than the convection. Note that 
two sets of curves are shown for the 35° cone. One (the solid lines) corre- 
sponds to entries at the 10 g acceleration limited boundary, and the other 
(dashed lines) to entries at higher altitudes. When the radiative heating is 
large, minimum total heating occurs at an altitude higher than that for the 
given acceleration and Reynolds number limits. This minimum is shown by the 
dashed curve labeled The ambient air densities at the flight altitudes 

are given in auxiliary scales along the abscissa. 

The total heating energy fractions for cone angles from 20° to 6o° are 
given in figure 10 , with the base radius, body density, and heat-shield prop- 
erties fixed. The minimum envelope of these curves gives the minimum value of 
heat input energy fraction that can be attained at each entry velocity for the 
conditions specified. The large angled cones depart from this minimum with 
increasing entry velocity and, as a result of their large radiative heat 
inputs, attain values of tj much larger than the minimum. The optimum cones 
have half -angles slightly less than 25 ° at V E = 3 V s and near 3 O 0 over a 
substantial part of the entry velocity range considered. 

Presenting the energy fraction as a function of cone angle, as in 
figure 11 , is useful for defining more precisely the optimum cone angle and 
the minimum energy fraction. In the case shown, which is typical of large 
limit Reynolds numbers, the minimum heat input within the given limits occurs 
at the intersection of Reynolds number -limited and acceleration-limited curves. 
Thus, the optimum cases attain the limiting values of both the acceleration 
and the Reynolds number at entry. On the dashed portions of the curves, one 
or the other of the limits is exceeded. Therefore, the minimum point in the 
Re z = 10- curve, for example, cannot be utilized. 

The definition of the optimum at higher and lower values of the limit 
Reynolds number is illustrated in figure 12. Curves are shown of r) as a 
function of 0 C for four sets of altitudes corresponding to limit Reynolds 
numbers from 10 5 to 10 s . At small cone angles, where radiative heating is 
negligible, t] varies on these curves inversely with (tan 0 C ) 3/2 according to 
equations ( 28 ) and (34). At larger cone angles, rj increases as the radiative 
heat input becomes important. Hence, each curve for a given Reynolds number 
limit exhibits a minimum. A fifth curve shows r) as a function of 0 C with 
altitude determined by the acceleration limit. In this case, tj is indepen- 
dent of Q c at small cone angles, where radiation is negligible (eqs. ( 28 ) 
and ( 33 ) )> rising as cone angle is increased to the region of important 
radiative heating. 
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Figure 10. - Variation of heat input energy fraction with entry velocity for cones of various 

angles on constant altitude trajectory. 
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Figure 11.- Variation of heating with cone Figure 12.- Effect of limit Reynolds number 


angle on constant altitude trajectory, 
Ve/V s =2.2. 


on optimum energy fraction and cone angle 
for V E = 3 V s . 
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Figure 13 . - Definition of minimum heat input 
and optimum cone angle for bodies of three 
base radii at Re^ = 10 7 . 


The values of minimum tj within 
the simultaneous limitations on accel- 
eration and Reynolds number are shown 
by the circled points for the four 
values of Re^. At Re^ = 10 8 and 10 7 , 
the optimum occurs in the region of 
small radiative heating, and is similar 
to that discussed in figure 11. At 
Re^ = 10 6 and 10 5 , the optimum has a 
significant radiative contribution, and 
occurs on the Reynolds number limited 
curve at an altitude above that for the 
acceleration limit. For example, at 
Re^ = 10 5 , the optimum occurs at a 
flight altitude greater than that for 
= 10 g. This is shown by the values 
of t] c , the dashed curves. (Note that 
tj c - 1 eq. (28).) The peak accel- 
eration at the optimum point for this 
Reynolds number limit is, in fact, only 
about 1.5 g. Hence, it is seen that 
the optimum point need not occur at 
accelerations as large as the accelera- 
tion limit. However, the Reynolds num- 
ber attains its limiting value in every 
case. Furthermore, lowering the limit 
Reynolds number below 10 7 can signifi- 
cantly increase the heat input and 
optimum cone angle, and also, by the 
arguments of an earlier section, the 
lift-drag ratio required at entry. The 
increase in r), however, is not as 
large as that which would occur with 
transition to turbulence in the absence 
of a Reynolds number limit. 

Curves like those in figure 11 are 
given in figures 13 for bodies of three 
base radii at entry velocities from 1.4 
to 3.0 V s , for a limit Reynolds number 
of 10 7 . These curves define the varia- 
tion of the minimum energy fraction, 
which increases with entry velocity, and 
the optimum cone angle, which decreases 
with increasing entry velocity. Compar- 
ison of figures 13 (a), (b), and (c) 
shows the base radius to be an important 
variable, and in figure l4, optimum 
energy fraction is given as a function 
of base radius for fixed Reynolds number 
and acceleration limits and body density. 
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Figure 14.- Dependence of optimum energy- 
fraction on base radius for constant 
altitude deceleration. 
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For small bodies, the optimum cone 
angle is small, the radiative heating 
is negligible, and the optimum energy 
fraction varies inversely with base 
radius. For large bodies, the optimum 
cone angle is increased, the radiative 
heating becomes large, and the optimum 
energy fraction increases with body 
size (dashed curves). An optimum base 
radius is defined. If, in place of the 
above constraints, the altitude is per- 
mitted to go above that dictated by the 
Reynolds number and acceleration limits 
to the altitude for minimum total heat 
input, the solid curves are obtained. 

On these, the heat input is essentially 
independent of base radius at large 
base radii. 

The high values of tj for the 
bodies of small base radius are a 
result of the acceleration limit 
imposed. Without this limit, the 
curves on the left of figure 13 ( c ) may 
be continued to their minimum points to 
obtain the low values of energy frac- 
tion shown in figure 15 with aim base 
radius. These values would be avail- 
able to unmanned probes. Note that the 
Reynolds number limit of 10 7 is still 
applied. 


Figure 15*- Energy fractions of a 1-meter base The body density also enters as a 

radius body without an acceleration limit, consideration. On acceleration-limited 

entries, tj c ~ 1/n/p^* Thus, the 

acceleration-limited curves of figure 13 ( c ) would be lowered a factor of 3 
the small cone angle range by a ninefold increase in body density. There is, 
therefore, an advantage in the use of dense vehicles in the small size range. 
With larger bodies (e.g., fig. 11), lowering the acceleration-limited curve 
allows the Reynolds number -limited curve to govern. The Reynolds number - 
limited cases are independent of body density. Hence, a small and limited 
reduction in r\ is available from body density increases at the larger base 
radii. 


The effect of the limit Reynolds number on optimum energy fraction is 
shown in figures 1 6 for large base radius entry bodies at several entry veloc- 
ities. The second ablation material of the analysis, Teflon, is introduced 
here for comparison. The characteristic effects of the Reynolds number limit, 
developed in connection with figure 12, are shown here explicitly. When the 
Reynolds number limit is high T] op ^ is relatively independent of Re^, but 
when the limit is small and the body base radius is large, T) op t is sensitive 
to Re^ and increases as Re^ decreases. The energy fractions remain below 
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(a) r h = 3.05 m (b) r b = 4.57 m 

Figure 1 6 .- Variation of optimum energy fraction with limit Reynolds number on constant 

altitude entries. 

0.017, however, for all cases shown, at Reynolds number limits as small as 10 6 
and entry velocities to 3 V s# Optimum cone angles for the conditions of fig- 
ures^ are shown in figures 17. The optimum cone angle decreases with increas- 
ing entry velocity and with increasing Reynolds number limit, and is relatively 
insensitive to ablation material. 

The low values of energy fraction calculated for Teflon (fig. l6) make 
that material appear attractive. These values are due, however, to the high 
blowing rates, that is, high mass loss rates, associated with the low heat of 
ablation of Teflon. Mass loss curves (figs. l8) , in fact, show the order of 
merit of the two materials and heavily favor the material with large heat of 
ablation. They also show the effect of changing the limit Reynolds number on 
the mass loss. The effect of the limit Reynolds number can perhaps be 




(a) = 3.05 m (b) r b = 4.57 m 

Figure 17*- Optimum cone angle variation with limit Reynolds number and entry velocity 

for constant altitude entry. 
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It is noted that if the asymptotic blockage factor a were taken to be 0 
instead of 0.1, the mass losses with Teflon at the higher entry velocities 
would be approximately halved, while those of quartz would be less affected 
(fig. l). Hence, the determination of the correct value of this constant is 
important. 

The other effect of low values of the limit Reynolds number referred to 
earlier is to increase the lift-drag ratio required. If the lift-drag ratio 
required exceeds that available, the entry becomes unattainable. The under- 
shoot altitude may then be said to have gone above the overshoot, so that the 
available entry corridor is negative. Conditions for which this may occur are 


shown in figures 19 , where L/d required at entry is plotted as a function of 



Figure 19* - Limitation on entry velocity imposed by joint consideration of Re^ and 

L/D available. 

Reynolds number limit for bodies of two (large) base radii. As long as the 
acceleration limit is attained at entry, the L/d is, according to equa- 
tion (30), independent of Reynolds number. This condition is represented by 
the horizontal sections of the curves. When the acceleration goes below the 
limiting value, the L/d needed rises rapidly, as on the curve sections at the 
left. The optimum cone angle also varies with Re^, as in figure 12, and along 
with it, the available L/D (appendix C). If the L/D required crosses the 
boundary marked maximum L/D available, it enters the unattainable region. 

This is seen to be important at the highest entry velocities considered, con- 
sistent with figure 3 , which shows that L/d becomes marginal at such entry 
velocities. 

Experimental observations of transition in very high enthalpy flow over 
ablating cones will be needed to determine a suitable value of the parameter 

Re z- 

Constant Reynolds number entries .- While the trajectory analysis showed 
that constant Reynolds number descents could be made with a 10 g acceleration 
limit, by designing to a descent rate of 25 m/s, the heating calculations were 
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not limited to such cases. Interest in the constant Reynolds number trajectory 
as a minimum heat input trajectory suggested its more general investigation. 


Curves similar in form to those presented for constant altitude entry were 
obtained. Thus, figure 20 shows the effect of entry velocity on the heat input 
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energy fraction for a given local 
Reynolds number (l0 7 ), cone angle, base 
radius, and ablation material. (Note 
that the energy fraction is independent 
of body density for the constant 
Reynolds number trajectories.) As with 
the constant altitude entries, where the 
radiative contribution is large, 

Reynolds numbers smaller than the speci- 
fied value can lower the total heat 
input. 


Figure 20.- Typical variation of the heat An ■ t 1]rp of > pn 

input energy fractions with entry velocity . An l nter estmg 1 eature . OX llgure 
on constant Reynolds number trajectory. is that the laminar convective energy 

fraction, rj , diminishes with increasing 

entry velocity. This will occur when a is small and K C V^ 
to 1. In the limit of o - 0, the integral in equation (28) 
pV - constant, 
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Figure 21 shows the minimum heat input energy fractions as a function of 
entry velocity for a series of cone angles at a given Reynolds number, and 
the envelope of minimum t] for these conditions is defined. Figure 22 shows 
the effect of body base radius and ablation shield properties on the envelope 
values. The ablation property effect is comparable to that on constant alti- 
tude entry, but the base radius effect is small. The convective heating con- 
tribution is, in fact, independent of r^, as may be seen from equation ( 28 ) 
where (pVr^) 1/2 may be taken outside the integral as a constant (see also 
eq. (B5) ) . Furthermore, the radiative contribution is only weakly dependent 
on r^, varying as r^ 0 * 2 according to equation (24) with p = 1.8. Hence, 
unlike the constant altitude entry, the constant Reynolds number entry shows 
heating energy fractions essentially independent of body size. 

The effect of the Reynolds number on the optimum energy fractions is 
shown in figure 23. The variation in the convective heat input with Re" 1/2 
is the principal part of this dependence, but changes with Reynolds number in 
the optimum cone angle and radiative input also contribute. The optimum cone 
angles are shown for four Reynolds numbers and two ablation materials in fig- 
ures 24 . There is some sensitivity to Reynolds number, while the effects of 
ablation material and body size are minor. In figure 25, mass loss fractions 
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Figure 24.- Optimum cone angles on constant Reynolds number trajectories. 


are presented as a function of entry 
velocity for optimum cones of the two 
ablation materials. 

Comparison of results for con stant 
altitude, constant Reynolds numb er , and 
ballistic entry . - The minimum energy 
fractions calculated for constant alti- 
tude entry, constant Reynolds number 
entry, and ballistic entry are compared 
in figure 26 for a common Reynolds num- 
ber limit, 10 Y , and a vaporizing quartz 
ablator. For consistency with the other 
results, the energy fractions on ballis- 
tic entry have been recomputed using the 
present heat-transfer relations. The 
constant Reynolds number descent has the 
lowest heat input, as was expected. For 
it and the ballistic entry, the energy fractions are insensitive to body size, # 
a consequence of the fact that neither is acceleration limited. For the con- 
stant altitude trajectories, the large sensitivity of energy fraction to body 
size noted earlier appears as a result of the imposed acceleration limit. The 
energy fractions of the acceleration-limited 1-m base radius bodies are, 
in fact, off the scale selected and not shown. The constant altitude entry of 
this size body without an acceleration limit, interestingly enough, falls 
within the band for ballistic entry. Hence, constant altitude entry without 
acceleration limit compares almost identically with ballistic entry in respect 
to total heat input. 


Teflon Vop. quartz 

Re € (million) = I 5 10 20 l 5 10 20 



Figure 25*- Mass loss fractions on constant 
Reynolds number trajectories. 


The optimum values in the best cases do not differ greatly for the three 
types of trajectories. It is noteworthy that grazing entries with accelera- 
tions within human tolerance can be made in large entry vehicles with heat 
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Figure 2 6 .- Comparison of optimum heat input energy fractions for three types of trajectories. 



input energy fractions and mass loss 
ratios comparable to the lowest possi- 
ble values obtained on either steep 
ballistic entry or the optimum, con- 
stant Reynolds number descents. 

Ablation of the vehicle nose . - The 
degree to which initially conical entry 
bodies become blunted by ablation dur- 
ing entry has been calculated for a 
number of constant altitude and ballis- 
tic entries by the method of appendix E. 
Noses of three materials. Teflon, vapor- 
izing quartz, and graphite, were consid- 
ered. One example is shown in figure 27, 
where nose radius to base radius ratio 
is plotted as a function of velocity 
for the quartz ablator on constant altitude entry of a 29 ° cone of an entry 
velocity of 2.6 V s . The rate of growth of nose radius with velocity, -dr n /dV, 
is nearly constant from entry down to V/V s = 1.8, after which it diminishes. 
The calculation is stopped at satellite velocity, at which point the rate of 
growth is small. Extrapolation to zero velocity gives a radius ratio less 
than 0.1, corresponding to a mass loss at the tip of the order of 1 percent 
of the total mass at entry. Furthermore, over 99 percent of the body surface, 
the bow shock wave is oblique during the period of intense radiative heating 


Figure 27*- Nose radius as a function of 
velocity for a quartz nose on constant 
altitude entry at Vg = 2.6 V s . 
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so that the radiative heat input is 
small in the presence of the small nose 
radius obtained. 


0; =lOg 
p b =320 kg/m 3 
Re,=l0 7 
O r b =3.05 m 
□ r b =4.57m 


Values of the radius ratio at 
satellite velocity are shown as func- 
tions of entry velocity in figure 28 
for constant altitude entries with the 
three tip materials and two body sizes. 
The final nose radii tend to become 
very large for Teflon ablators, and 
moderately large with vaporizing quartz 
at the higher entry velocities, but 
remain less than 0.1 r^ for graphite. 

1.8 2.0 2.2 2.4 2.6 2.8 3.0 

For ballistic entry comparable 

Figure 28.- Nose radius at satellite velocity results not calculated previously are 
as a function of entry velocity for ini- shown in figures 29 and 30 for bodies 

tially sharp cones of three materials on with a base radius of 1 m. The ablation 

constant altitude entry. Q f a Teflon nose for Ve/V s = 2 . 0 , 

figure 29 , has a history much like that in figure 27. Figure 30, which 
extends to higher entry velocities than figure 28, shows that on ballistic 
entry as well as on grazing entry, nose ablation can be severe or moderate, 
depending on the nose material. 

While these estimates should not be regarded as definitive (in view of 
the. assumptions, appendix E), they suggest that optimum conical entry bodies 
can substantially retain the advantages of pointed cones with respect to 




0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 

_v_ 

v s 


Figure 29- - Nose radius as a function of 
velocity for a Teflon nose on ballistic 
entry at an entry velocity of 2 V s * 
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Figure 30. - Nose radius at satellite speed as 
a function of entry velocity for initially 
sharp cones on ballistic entry. 
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radiative heat input, 6 and can have acceptable mass losses at the tip, 
provided high performance ablation materials, such as the assumed vaporizing 
quartz and graphite, can be successfully applied. 


REVIEW OF AS STJMPTI ON S 


The present analysis was performed to demonstrate the existence of 
optimum conditions, to define the magnitude of the optimum heat input, and 
cone angles, etc., and to identify the effects of the principal variables in 
parametric fashion. By its nature, such an analysis of a complex problem must 
be approximate, but the simplifying assumptions made to permit analysis 
must not invalidate the results to be obtained. To emphasize this, it is use- 
ful to review and comment on the assumptions made, which are as follows: 

1. The angle between the flight path and the horizontal is small. 

Comment: Satisfied exactly for constant altitude entry, and very well satis- 

fied on high-speed portions of constant Reynolds number entries. 

2 . Changes in mass, cone angle, and base radius of the vehicle are small 

during entry and may be neglected. Comment: The cone angle and base radius 

changes are indeed small for cases considered. The mass loss integrated over 
the trajectory is restricted to 0.2 of the mass at entry or less in cases 
presented. 

3. The drag coefficient is constant during entry. Comment: Should be 

valid in hypersonic speed range for bodies of the class studied. 

4. Laminar boundary layer can be maintained by restricting the local 

Reynolds number below a suitable limit. Comment: Assumption is consistent 

with knowledge of boundary-layer transition processes for nonablating surfaces 
at lower speeds. Its validity at the speeds and ablation conditions of this 
analysis requires experimental study. Considerations of roughness of the 
ablating surface, uniformity of vapor emission with time and over the surface, 
etc., may enter. 

5* The ratio of the Reynolds number evaluated at boundary-layer edge 
properties to that evaluated at free-stream properties is a constant for any 
one angle. Comment: Approximately true, and, as shown in appendix D, the 

constant chosen will yield a local Reynolds number equal to or less than the 
selected limit. 

6. The density profile of the earth’s atmosphere has an exponential 
dependence on altitude. Comment: Not a critical assumption for this analysis. 


6 A favorable effect of tip blunting is to reduce the convective heating 
rate. The reduction is sizable when the shear layer of the body nose becomes 
larger than the boundary-layer thickness. This effect has (conservatively) 
been neglected in the present analysis. 
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7. The heat-transfer coefficients may he represented by equations given 

in the text. Comment: Much of the speed range considered is beyond the range 

of present day laboratory tests. The equations employed are based on pres- 
ently available theoretical and experimental information, and use best esti- 
mates of high-temperature gas properties. For conical flow, less extrapolation 
of radiative and convective transport properties is required than for stagna- 
tion point flow, since the enthalpy is well below stagnation enthalpy. For 
example, the radiative heating on optimum cones at V;g/V s = 3 is within the 
range of present day experiments on blunt bodies, and should not be in doubt 
by more than a factor of 2. With larger cone angles and blunt bodies, larger 
uncertainty exists in the radiative transfer. For example, the vacuum ultra- 
violet contribution, recently discovered to be important, has not been included 
in the formulas presented, and would increase the heat input for U >1.4 V s . 
Convective heating is presently believed to be subject to less uncertainty 
than the radiative heating in this high-speed range. 

8. Flow along streamlines outside the boundary layer may be regarded as 

isoenergetic, that is, flow energy losses by radiation may be neglected. Com- 
ment: The heat-transfer coefficient is defined in terms of the available flow 

energy. Small values of the radiative heat -transfer coefficient therefore 
signify essentially constant energy flow and, from the results presented, it 
is clear that the optimum bodies satisfy the assumption. Nonoptimum bodies, 
however, may have large values of Cji e , and deviate from isoenergetic flow. 
Several cases were analyzed allowing for energy loss along streamlines and one 
is shown in figure 31. The energy input is affected by only a few percent for 
this 60° half-angle cone. For smaller angle cones, the effect is even less. 


9* Radiative reabsorption in the shock layer may be neglected; that is, 
the shock layer is assumed to be optically thin. Comment: If radiative inten- 

sity at all wavelengths is small compared to that of a black body of the shock- 
layer temperature, the assumption is justified. This appears to be the case 
for the optimum and near-optimum cones analyzed. For large angled cones with 
shock-layer temperatures high enough to produce significant ionization, 



Figure 31. - The effect of a nonadiabatic shock 
layer on the heating of an effectively 
sharp cone of 0 = 6o° on constant 
altitude entry. 


U >1.4 V s , the vacuum ultraviolet 
region of the spectrum becomes impor- 
tant (refs. l6, 17 ) and reabsorption of 
radiation may be an important consider- 
ation. Other atmospheric gases, notably 
those generating CN as a radiating 
species, can also encounter important 
reabsorption at particular wavelengths. 
With air, no such predominant radiators 
are present. 

10. The heat blocking and ablative 
behavior of the heat shield may be sim- 
plified, as represented by equations ( 9 ) 
and (l8). Comment: This assumption has 

four parts. The form assumed for the 
heat blocking factor ^ has been dis- 
cussed in the Analysis and Presentation 
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of Results. Although the magnitude of the heat input is affected by the value 
assigned to the asymptote cr, the variation of heat input with cone angle is 
unaffected, and the existence and nature of the optimum do not depend on the 
value chosen. Second, chemical reactions at the surface and in the boundary 
layer have been neglected. For the quartz material, this is appropriate, but 
for graphite and Teflon, some chemical reaction must be expected. However, at 
the very high speeds considered, reaction products tend to be dissociated at 
boundary-layer temperatures, and chemical reactions may not be of overriding 
importance. Third, heat reradiated by the ablation shield is neglected. Com- 
parison of the rate at which heat can be reradiated with the heat input rate shows 
that the reradiation is small compared to the heat input rate for surface tem- 
peratures below about 3500° K. This is near the maximum permissible surface 
temperature for any material, so the assumption may be allowed. Fourth, the 
production of ablation vapor is assumed to be a result of convective heating 
only. For cones of near -optimum angle, this is a good approximation. For 
cones of larger angle experiencing large radiative heat input, the transpira- 
tion rates would be higher than calculated, with a moderate reduction in con- 
vective heat input. It is clear that the treatment of the ablation shield is 
approximate, and that some basic understanding of the high velocity behavior is 
still lacking, but it should be adequate to indicate the mass loss and the heat 
blocked to a first approximation. 

11. The heating of partial cones, modified to generate lift, may be 

approximated by analysis of the complete cones from which they are derived. 
Comment: Intuitively, one can see that this will be the case to a useful order 

of approximation. In detail, a number of corrections or modifications to the 
analysis would be required to treat the modified cone. For example, in the 
case of oblique base plane modification described in appendix C, if the mean 
base radius of the upper and lower meridians is regarded as the base radius of 
the analysis, the convective heat input is reduced on the upper meridian and 
increased on the lower meridian in approximately compensating amounts. The 
same is true for radiative input. The Reynolds number limit must be applied 

to the longest (lower) meridian, which decreases the stream density permitted, 
and is perhaps the main effect. Certainly, the magnitude of the over -all 
result will not be greatly affected. 

12. The heating increment on tangential entry into the atmosphere from 

the point of entry to the point of joining the selected constant altitude or 
constant Reynolds number trajectory is equivalent to that obtained by entering 
the design trajectory at full entry velocity. Comment: For vehicles of the 

size, density, etc. , considered herein, velocity losses of the order of a few 
percent of Vg will, in actual entries, occur before design altitude is 
reached. Since the Reynolds numbers are, in this phase, lower than the design 
value, convective heat input will be increased (eq. ( 28 )), while radiative 
input will be diminished (eq. (22)). The altitude range of the principal 
velocity loss is within one atmospheric scale height of the design altitude, 
so that the convective heat input may be shown to be increased by a factor 
smaller than e 1/2 . If, in fact, a factor of 2 is allowed on the values of 

Cji for this interval, the increment in r\ is typically 0.0002. The 
approximation is therefore satisfactory. 
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CONCLUDING REMARKS 


The results presented indicate that heat input energy fractions on hyper- 
bolic grazing entries, -with acceleration limited to within human tolerance, can 
be limited to about the same levels as on ballistic entries which have no 
acceleration limit. Under the best conditions, corresponding to moderately 
high boundary-layer transition Reynolds numbers, entry may occur at velocities 
up to 3 V s with less than 1 percent of the initial kinetic energy finding its 
way into the body as heat. Corresponding mass losses of the heat shield may 
be less than 20 percent of the total mass at entry with shield materials of 
high heat of ablation. Conical bodies away from the optimum, of which blunt 
bodies may be considered one extreme, will experience heat inputs many times 
greater than the optimum. 

The optimum bodies are generally swept back to such an extent that the 
radiative heat input is small compared to the convective. An exception to this 
arises if the local Reynolds number must be limited below one million to main- 
tain laminar boundary layer, in which case the radiative input may comprise the 
order of one -third the total heat input at optimum. The total heat input is 
substantially greater in these cases than in the case of higher limit Reynolds 
numbers. 

The determination of the level of Reynolds number for maintaining a 
laminar boundary layer is an important problem requiring experimental solution. 
Not only are the heat input and heat -shield mass loss increased for low values 
of the Reynolds number limit, but the acceleration may be reduced well below 
its limiting value with the consequence that a high lift-drag ratio is 
required. At values of entry velocity approaching 3 V s , the L/D required 
tends, in any case, toward higher values than can be aerodynamically generated 
by modification of the optimum cones. Lowering the total acceleration to 
obtain a low Reynolds number will, at such velocities, usually insure that 
insufficient aerodynamic lift is available. 

For entries into the earth’s atmosphere at or below 10 g total accelera- 
tion and with levels of L/D which can be generated with cones of optimum 
half -angle, the maximum entry velocity will be limited to about 3 V s . Higher 
velocity entry might be accommodated by use of multiple pass entries or by use 
of rocket thrust to supply lift or drag force. 

The trajectories selected for analysis tend to give near the minimum 
possible convective heat input by staying near the highest permissible Reynolds 
number. Trajectories on which the Reynolds number is allowed to go far below 
this limit will increase the heat input, as will those that go above it and 
incur a turbulent boundary layer. The constant Reynolds number descent obtains 
the minimum convective heat input exactly, for a given limit Reynolds number, 
cone angle, entry velocity, and ablation shield. If on this trajectory the 
cone angle is selected to make the total heat input a minimum, and if the 
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radiative input is found to be small, then the entry may be identified as a 
minimum heat input case with respect to both trajectory and cone angle. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 7 > 1966 
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APPENDIX A 


THE LAMINAR CONVECTIVE HEAT -TRANSFER COEFFICIENT 


The problem of laminar convective heat transfer to nonablating cones at 
very high flight velocities has been investigated theoretically by G. Chapman 
(ref. 10). Chapman finds that his many solutions for a variety of flight 
velocities from 4.5 to 25 km/s, cone angles from 15° to 60° , and wall temper- 
atures from 1000° K to 4200° K can be correlated by equations of the form 


Nu = F 

HZ . r 0.36 

^Re w V 


where 


Nu 


Re w 


1 c Pw X 

h t - K K 

Pw u e x 

M-w 


(Al) 


(A2) 

(A3) 


and q is the local heating rate at a distance x from the cone vertex. 1 * 
The coefficient F is weakly dependent on wall temperature and surface pres- 
sure (see fig. 3, ref. 10), and takes the following values when V is in 



Tw, °K 
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8.45 
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10.11 
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Equation (Al) may be put in the form of the heat-transfer coefficient used 
herein if the heat -transfer rates over the surface are integrated and the 
integral is equated to the present defining equation of total heating rate at 
zero ablation: 


Htt pV 3 A _ 4 F ht hv n/'p'W j, w u £ a 

Hc o 2 3 y°« 36 ^r^ sin Q c Rr w 


(A4) 


1 Some departures from the correlation occur when wall enthalpy approaches 
recovery enthalpy and when the boundary layer becomes highly ionized. 

Chapman (ref. 10) discusses the extent of these deviations and the conditions 
under which they occur. It is sufficient here to note that they do not 

importantly affect the optimum conditions of the present analysis. 
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so 


For adiabatic floy h^ = h + 
h^ - h^. may be replaced by 



V 2 / 2, and furthermore (h - h^j « V 2 /2, 
V 2 /2. Hence, 

/ £2. n HE 

4 F v P V 

3 y °*86 s -j_n 0 C Pr w 


(A5) 


In addition, u e /V may be replaced by cos G c to a good order of approxima- 
tion, and 


Py _ T 
P “ P Z y T y 


V 2 sin 2 0 C 
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+ 1 
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yhere the yall pressure ratio is converted to the terms In the bracket by use 
of Neytonian theory, and V 2 sin 2 0 c /RT » 1 for the velocities and cone 
angles considered. Equation (A5) becomes 
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yhich is of the form of equation (l6) yith 

, r * 


c = - - 
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Z W RT- 


Pr. 


'y nJ -y 
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(A7) 


From this, it may be seen that C c depends on surface temperature and pres- 
sure, through F and through the air properties at the yall . In neither case 
is the dependence a strong one, hoyever. For example, ^ ]± M / T v increases by 
only 27 percent as T v is lovered from 3000° K to 1000° K, yhile Z w remains 
very close to 1.0. The conditions chosen for evaluating C c were T w = 3000°K 
and p w in the range from 0.1 to 1.0 atmosphere. With these selections, 


/ V \°’ 14 

C c = I.l6xl0" 4 l 


(A8) 


Comparison of equations (l6) and (17) yith the laminar convective heating 
equations of reference 5 shoys that in the loyer speed region of that refer- 
ence, V < 13,000 m/s, the given values of c H Co agree to yithin about 10 per- 
cent. The higher speed equation of reference 5 "was based on extrapolation of 
a loyer speed theory, and underestimates the present results by factors ranging 
from 1.1 at 13,000 m/s to about 4.0 at V = 3°; 000 m/s. This fact should be 
noted In comparing numerical results of reference 5 yith the present results. 
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APPENDIX B 


INTEGRATION OF CONVECTIVE HEATING RELATION FOR 
CONSTANT REYNOLDS NUMBER TRAJECTORIES 


For the constant Reynolds number trajectory the laminar convective heat- 
transfer relation, equation (28), may he integrated in closed form for j = l/6. 

„ . c °^ ln . 28 ° . r Y *( (BD 

VjfJZ sin 2 0 C ^o ' 1 + K C V 2 / / JL 

^ Po V 

where, for constant Reynolds number, pV = K is constant. The integral may 
be expressed as a sum of two integrals 
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The first integral can be evaluated immediately, while the second requires the 


change of variable Y 


= W2/3 


, which results in 
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The solution to this integral can be found in reference l8 and is 
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After evaluating this expression at the limits of integration and returning to 
the original variable V, we find 
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APPENDIX C 


LIFT -DRAG RATIOS OF MODIFIED CONES 


Three -ways of generating lift with cones are known or have been proposed: 
(l) controlling attitude of the unmodified cone; (2) removing all or part of 
the upper conical surface above a given plane surface; and (3) using an oblique 
base plane. The latter two alternatives are sketched in figure 32. 


The first alternative is the least 
attractive for hyperbolic entry, since 
it creates crossflow conditions in the 
boundary layer that are very damaging 
to laminar stability (see, e.g., 
refs. 19 and 20). Furthermore, since 
the effective angle presented by the 
lower surface to the stream is G c + a t , 
and since L/d and a-^ must diminish 
with velocity, the effective cone angle 
and drag coefficient diminish during 
entry increasing the total heat input. 

At any speed, Cp is less than its 
allowable maximum which is dictated by 
the angle presented on the lower sur- 
face only. Hence, this alternative does 
not recommend itself compared to the 
other two. 

Figure 32.- Lift-drag ratios of modified 

cones at zero angle of attack. ^ ie second alternative, proposed 

in reference 21, provides the maximum 

lift-drag ratio in the limiting case of the half cone. That ratio is shown 
as a function of cone angle in figure 3^* Any smaller value of L/D down to 
zero may be selected by canting the plane of the upper surface at an angle 
facing the stream. In the absence of skin friction and base drag, this modi- 
fication gives infinite L/d at 8 C = 0. Cases of present interest, for 
example, 0 C = 30°, have maximum l/d up to 1.1. Crossflow affects this con- 
figuration along the intersection of the plane surface with the conical sur- 
face, where a pressure jump exists laterally (in the case of a rounded 
intersection, distributed to give a lateral pressure gradient). Turbulent 
boundary layer might be expected to appear first along this region and to 
spread over the adjoining surfaces at appropriate Reynolds numbers. The con- 
vective heating of the upper surface imposes a penalty not considered quanti- 
tatively in the analysis. The estimation of the heating of the upper surface 
poses some theoretical problems in the presence of a blunted tip and the 
crossflow phenomenon described above. 

The third alternative, proposed in reference 22, was investigated by 
means of Newtonian theory in reference 23 ^ and the results for lift -drag ratio 
are also shown in figure 32. The lift-drag ratio of the cone with an oblique 
base plane is the same as for a flat plate inclined at the base angle to the 
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stream, that is, L/d = cot 6 . As 5-^0, L/d becomes very large. However, 
for bodies enclosing a real physical volume, 5 must be greater than 0 C , so 
that for e c = 30 °, a base angle of 42° would comprise a very thin and doubt- 
fully practical case, which would have the same L/d as the 30 ° half cone, 
namely, 1.1- Bodies of this cone angle with a useful volume arrangement might, 
in fact, be restricted to lift-drag ratios of 0.8 or less. This modification 
should, however, be free of crossflow effects, except those occurring behind 
the vehicle base, and may present the best possibility for retaining laminar 
flow over the conical heat shield. It is also the closest of the three alter- 
natives to the assumptions of the heating analysis, which treats a pure cone 
at zero angle of attack. 

With alternatives (2) and ( 3 ), the method of reducing the component of 
l/d in the vertical plane as velocity diminishes along the trajectory is to 
roll the vehicle until L/D times the cosine of the roll angle equals the 
L/D required. To prevent lateral curvature of the flight path, a roll oscil- 
lation symmetrical about the position for vertical lift vector is followed, 
the amplitude determining the vertically effective component of L/D. The 
angles of pitch and sideslip of the vehicle are (ideally) zero throughout 
entry, simplifying the theoretical treatment of heating. 
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APPENDIX D 


DEVELOPMENT OF TRAJECTORY RELATIONS 


In this appendix, the first section develops the equations for 
acceleration-limited, Reynolds number -limited constant altitude flight; the 
second develops the equations for constant Reynolds number descent; and the 
third defines the relationship between local and free-stream Reynolds numbers 
on cones in hypervelocity flight. 


EQUATIONS OF CONSTANT ALTITUDE DECELERATION 


The lift force required to maintain circular flight above satellite 
velocity is given by 


L = [(V/V s ) 2 - 1]W 

from which it follows that the lift-to-drag ratio is 


L 

D 
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(Dl) 
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For zero dynamic pressure q , L/D is required to be infinite^ since the need 

00 / 

for the lift force remains, while the drag has gone to zero (note asymptotic 
behavior of the curves in figure 3) • To fully display the effect of velocity 
on L/D, equation (D2) is written 


L = (V/V R ) g - 1 ^/CpA , 

D (V/V s ) 2 (1/2)PV S 2 

where the factor on the right is a constant for a given vehicle and altitude, 
and may be identified as the L/D required at infinite velocity with this 
vehicle and altitude. 


The altitude is selected to maintain specified limits on the total 
acceleration and Reynolds number. The axial acceleration may be written as 


a x (l/2)pV s 2 / v \ 2 
g W/CpA 


(d4) 


while the component radial to the earth is the centripetal acceleration 
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so the total acceleration is 
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By equation (D3), this may he written entirely in terms of V and L/d. 
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which is solved for L/d to yield equation ( 30 ) of the text. The maximum 
value of a/g occurs when V = Vg and, in acceleration limited cases, is set 
equal to the limit value of acceleration, a^/g. The L/D varies after entry 
with velocity according to equation (D 3 )j while the acceleration decays accord- 
ing to equation (D 6 ). The stream density at which the peak acceleration 
attains the limit value may be obtained from equation (d 6 ). 
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If, on the other hand, the altitude is determined by the Reynolds number 
limit, then at entry onto the constant altitude trajectory, when Reynolds 
number is maximum. 
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where the subscript e identifies fluid properties at the boundary -layer edge 
and S is the slant length of the cone. As will be shown, to a good order of 
approximation, in the following equation 


pVr b /|i tan 9 C 


(DIO) 


k = k(u) only. The parameter k(u) has an approximately constant value of 2 
for U < V s and diminishes to 1.3 at U > 1.3 V • Through equation (DIO), 
equation (D9) may be written in terms of free -stream variables as 
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and the free-stream density determined by the Reynolds number limit is 
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P = Re Z 


tan 6 C \i 
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Thus, with the density for the acceleration, limit defined by equation (d 8 ) and 
for the Reynolds number limit by equation (D12) , the value selected to observe 
both limits is the smaller of the two. 


EQUATIONS OF CONSTANT REYNOLDS NUMBER DESCENT 

The condition that p £ u_S/p e remain constant along the trajectory 
requires that pVr^/n be approximately constant also. If we assume that 
variations in the stream viscosity with altitude may be neglected, and that 
iq, is not appreciably diminished by ablation in flight, then by equation (DIO) 


pV = Re g ~ ba ^ ~ = constant = K (D13) 

and equation ( 5 ) which, for small path angle 7 , governs deceleration along 
the trajectory becomes 


so that 


1 dV K 

V dt "m2 


V 

V E 


C P A k 

m 2 


e 


Similarly, in terms of x as independent variable, 


and 


d V = 2^ K 

dx m 2 


v = v E 


C P A k 

m 2 


X 


(D14) 


(D15) 


(di6) 


(D17) 


The rate of descent through the atmosphere implicit in equations (D13), 
(D15), and (D17) may be obtained as follows. Assume, following reference 1, 
that 



(D18) 
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By (D13) and (Dl 8 ) , 


p e" Py V = K 
K o 


Eliminating V between (D15) and (D19) gives 


(D19) 


or 


K 

P C V E 


n C D A K 

2 * 


= 1 


Py + 


CpA K 
m 2 


Zn 


*o V E 

K 


Zn 


Po V E 

pV 


(D20) 


(D21) 


Equation (D21) shows that the rate of descent is steady, as can be made clear 
by taking the derivative to obtain 


V = d y c d a k 

y “ dt m 2 p 

The flight-path angle is obtained as a function of the velocity from 


(D22) 



sin 7 


V y CpA K 

V m 2pV 


(D 23 ) 


and may also be obtained as a function of 
t or x by use of (D15) or (D 17 ). 

The lift -to -drag ratio required to 
follow this path and the accelerations 
experienced may be derived as follows. 

The motion is not circular about the 
center of the planet, and to find L/D, 
we must consider the actual flight-path 
curvature. The notation employed is 
indicated in sketch (b). About the 
instantaneous center of curvature 


R d/' = dx (D24) 


The angle 7 1 is measured relative to 
the reference direction C-l, while the 
local flight -path angle 7 is referred to 
the local tangent to the planet T s surface. 
From the sketch 


(D25) 


Sketch (b) 


d 7 1 = d 7 + dg 



Since d£ is a differential, we may write 

Rp d| = (dx)cos 7 

which is combined with (D24) and (D25) and solved for R to obtain 


(d//dx)+ [(cos 7 )/R p ] 
By Newton 1 s second law of motion 


tt t tt 2 f dy cos 7 

W cos 7 + L = mV 2 


from which ^ noting that gRp = V s 2 , 


L f V 


w vv 


R p tx + cos y ) ■ cos y 


By equations (D23) and (D17) 


sin 7 = 


CdA K/20 


m V E - (CpA/m) (K/2)x 


and by differentiation of (D 30 ) 


dy 

dx “ 


P sin 7 tan y 


so (B29) may be written for the constant Reynolds number descent 


— = J s ^ n 7 tan 7 + cos 7 j “ cos 7 


By the definition of drag coefficient 


0 , hi i pV 2 

W W 2 p 


¥ 2 V 


and L/D is obtained by dividing (D32) by (D33) 


L [(V/V s ) 2 - l]cos 7 + (V/V S ) 2 0R P sin 7 tan 7 


D 


(C D A/W)(KV S /2)(V/V S ) 


(D2 6) 


(D27) 


(D28) 


(D29) 


(D30) 


(D31) 


(D32) 


(D33) 


(D34) 
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For small 7 , 


cos 7 = 1 , 


sin 7 = tan 7 = Vy/V, 


and we may write 


L 

D 


_W 2_ Vs 

CdA KV s V 



- 1 + PRp 



(D35) 


With Vy = 0 , this reduces to the expression for constant altitude flight, 
equation (D2). With nonzero Vy, a further reduction is accomplished hy 
relating the factor outside the bracket to Vy/V through equation (D23)« 


L = 1 - (V S /V ) 2 | Jy 

D PRpVy/V + V 


(D36) 


Thus, the L/D required is obtained as a function of velocity and descent 
rate, while the descent rate has a known dependence on the vehicle parameters 
and the constant Reynolds number selected. 

The axial acceleration is given by equation (d 4). The transverse acceler- 
ation for small 7 is, by (D27), (D28) , and (D3l ) , 



(D37) 


which may be compared with equation (D 5 ) • 


The total acceleration is given by 


a 

g 



KVa 


L 2(V/V S )(W/C D A)- 



(D38) 


RELATIONSHIP BETWEEN LOCAL AND FREE-STREAM REYNOLDS NUMBER 


The ratio of the local Reynolds number based on edge properties and slant 
length to the free-stream Reynolds number, referred to base radius, may be 
written 

P € u e ^ S 

Re PVH e r b 

For a thin hypersonic shock layer 

u e 

— = cos 0 C 



and from the cone geometry, S/r^ = l/sin 0 c . Furthermore, 

Pe . p(V sin 0 C ) 2 
P e " RT e Z 6 - RT € Z e 
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2.4 


so equation (D39) "becomes 



Figure 33*- Ratio of local to free-streaxn 

Reynolds number for cones in hypervelocity 
flight. 


Re G (v sin 0 C ) p, cos 0 C 
ReT = RT e Z e jj^“ sin~0 c 


(d4o) 


In figure 33> bb e ratio (Re e /Re)tan 6 C 
is plotted as a function of normal 
velocity component, V sin 0 C , over the 
range from 5 to 12 km/s for altitudes 
ranging from 42*7 bo 6l km. The thermo- 
dynamic properties were evaluated with 
the aid of reference 24, and refer- 
ence 25 was used for the transport 
properties. 

The flight conditions shown in 
figure 33 are representative of those 
for near-optimum bodies. If Reynolds 
number ratio is written as 


Re 

— tan 0 C = k 


(d4i) 


then, k = 2 for normal velocities from 5 to 8 km/s and decreases to about 1.3 
at U = 10 km/s. In order to be conservative, and for convenience, k = 2 
has been used throughout the analysis. Use of this value of k at the higher 
values of U has the effect that for a given value of Re, Re € is less than 
the limit Reynolds number which is being imposed. 
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APPENDIX E 


ANALYSIS OF NOSE ABLATION 


The heating rate per unit area at the stagnation point of a blunted nose 
may be written 


<lH s t 

dt 


- °H st \ ev 3 


(El) 


where ^H s t ^l^e^lonless local stagnation-point heat-transfer coeffi- 

cient consisting of a laminar convective and a radiative component. The con- 
vective component may be written, based on theoretical results given in 
reference 26, as 


CH s t c = ( 3 • ^5x10 _4 //pr^) 


(E 2 ) 


where the constant 3.65XIO" 4 has the dimensions (kg/m 2 ) 1 '" 2 , and \|r represents 
the reduction in heating due to transpiration and is assumed to be given by 
equation (l8) • 

The radiative component may be written, assuming the shock layer is 
optically and physically thin. 


CH s t e = C e SV a ' 3 ill (E 3 ) 

p o 

where C e is given in the text and 5 is the shock-wave standoff distance 
at the stagnation point, given for spherical noses (ref. 27) by 

8 = O.78 r n - 2 - (E 4 ) 

Here, p is the gas density behind the shock wave, and is, for equilibrium 
shock layers, a function of p and V only. The heat-transfer coefficient 


C H s t = C H s t c + °Hst e 

is thus a function of ambient density, velocity, and nose radius. The rate 
of mass loss at the stagnation point may be written 


m lH s t/dt 
A = 5 


(e6) 
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and the rate of surface regression is 


so that, by (e 6) and (El), 



x n = 


^H st 2 


~ pV" 




(E7) 


(E8) 


The following assumptions are now made: (l) The nose is originally a 

sphere segment of small finite radius and remains spherical during ablation. 

(2) Deviations of the heat transfer from that given by the above relations, 
due to deviations from chemical and thermodynamic equilibrium, are small. 

(3) Energy release and surface erosion due to chemical reactions of the abla- 
tion materials with atmospheric gases may be neglected. (4) No material is 
removed by spallation. 


Under the assumption of spherical blunting, the rate of increase in nose 
radius may be related to the rate of surface regression at the stagnation 
point 


n 


x. 


■n 1 


sin 0 C 
- sin 0 C 


(E9) 


so that 


. _ C H s t g sin 0 C 

£p b 1 - sin e c 


r n = 


(E10) 


Given the trajectory relations, which define the variations of p and V with 
time, and the heat-transfer coefficient, equations (E2) through (E5), equa- 
tion (E1Q) may be numerically integrated to obtain the nose radius as a func- 
tion of time during entry. 
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